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Abstract 

We prove central and non-central limit theorems for the Hermite variations of the 
anisotropic fractional Brownian sheet W a ^ with Hurst parameter (a,/3) G (0, l) 2 . 
When < a < 1 — ^ or < /3 < 1 — ^ a central limit theorem holds for the 
renormalized Hermite variations of order q > 2, while for 1 — i < a, 8 < 1 we prove 
that these variations satisfy a non-central limit theorem. In fact, they converge to a 
random variable which is the value of a two-parameter Hermite process at time (1, 1). 
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1 Introduction 

In recent years a lot of attention has been given to the study of the (weighted) power 
variations for stochastic processes. Let us recall the case of the fractional Brownian motion 
(fBm). Consider B H := (-B^) te j ,i] a fBm with Hurst parameter H in (0,1). Recall that 
B H is a centered Gaussian process with covariance Rii(t, s) = E(B^ Bf) = \{t 2H + s 2H — 
\t — s\ 2H ) for every s,t € [0,1]. It can also be defined as the only self-similar Gaussian 
process with stationary increments. The weighted Hermite variations of order q > 1 of B 
are defined as 

N 

V N ■= ^2 f( B (i-l)/N) H q( NH \ B i/N ~ B fi-l)/N))-> 
i=l 
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where H q denotes the Hermite polynomial of order q (see Section I2.2|) and / is a real- valued 
deterministic function regular enough. Take for example q = 2. Since the second-order Her- 
mite polynomial is H 2 (x) = ^-g— 3 the latter quantity is equal to \ Yl!i=§ f(.B(i-i)/n)(N 2S l-^i/N' 
B^_iy N \ 2 — 1)- The asymptotic behavior of these variations plays an important role in es- 
timating the parameter of the fractional Brownian motion or of other self-similar processes 
(see e.g. [7] or [S])- Weighted Hermite variations are also crucial in the study of numerical 
schemes for stochastic differential equations driven by a fBm (see [12J ) . A full understand- 
ing of V n is given in j3j El EJ [21], when / = 1 and in p~U [T7] for quite general functions /. 
Let us recall the main results for the case / = 1: 



If < H < 1 - l/(2q) then, N^Vn Af(0, C H ), 

n— >oo 

If H = 1 - 1/(2?) then, (log^)^)- 1 / 2 ^ ^ M(0, C H ), 

n— >oo 

If 1 - l/(2g) < H < 1 then, N l - 2H V N ^ Hermite r.v. 



Here Ch is an explicit positive constant. A Hermite random variable is the value at time 
1 of a Hermite process, which is a non-Gaussian self-similar process with stationary incre- 
ments living in the qth. Wiener chaos (see e.g. [14]). 

In this paper we use Malliavin calculus and multiple stochastic integrals to study the asymp- 
totic behavior of the non- weighted (i.e. / = 1) Hermite variations, where the fBm is replaced 
by a fractional Brownian sheet (fBs), which is a centered Gaussian process (I^("'f))(s,t)G[o,i] 2 
whose covariance function is the product of the one of a fBm of parameter a in one direc- 
tion and of the covariance of a fBm of parameter f3 in the other component. We define the 
Hermite variations based on the rectangular increments of W a, P by, for every N,M > 1, 

N-1M-1 , / xx 

V ■= E E H » r ^ kSf i±i - Wfm Wf ± . 

.„ . „ V V N ' M N' M N ' M N'MJ J 

1=0 j=0 

In some sense the fBs is the tensorization of two orthogonal fBm. In the view of results 
from the one-dimensional case mentioned above it would be natural to expect the limit of 
the correctly renormalized non- weighted Hermite variations to be the " tensorization" of the 
limits appearing in the one-parameter case. This actually is true in several cases: when 
a, /3 < 1 — the limit in distribution of the renormalized sequence Vn, is, as expected, a 
Gaussian random variable and when a, ft > 1 — ^ the limit (in L 2 actually) of Vn, is the 
value at time (1, 1) of a two-parameter Hermite process (which will be introduced later in 
the paper) . The most interesting and quite unexpected case is when one Hurst parameter is 
less and the other one is strictly bigger than the critical value 1 — . It turns out that in this 
situation the limit in distribution of the renormalized Hermite variations is still Gaussian. 
We prove our central limit theorems using Malliavin calculus and the so-called Stein's 
method on Wiener chaos introduced by Nourdin and Peccati in [15j . Using these results, it 
is actually possible to measure the distance between the law of an arbitrary random variable 
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F (differentiable in the sense of the Malliavin calculus) and the standard normal law. This 
distance can be bounded by a quantity which involves the Malliavin derivative of F. Using 
these tools and analyzing the Malliavin derivatives of Vn m (which is an element of the q 
th Wiener chaos gener ated by the ffis W a > p ) we are able to derive a Berry-Esseen bound in 
our central limit theorem. 



We proceed as follows. In Section[2]we define the Hermite variations of a ffis and give 
the basic tools of Malliavin calculus for the fractional Brownian sheet needed throughout 
the paper. The central case is presented in Section EJ whereas Section [4] is devoted to the 
non-central case. The Appendix contains some auxiliary technical lemmas. 



2 Preliminaries 

2.1 The fractional Brownian sheet 

Several extensions of the fractional Brownian motion have been proposed in the literature 
as for example the fractional Brownian field ([HJH]), the Levy's fractional Brownian field 
([6]) and the anisotropic fractional Brownian sheet ( [TOj [2] ) , which we consider in this paper. 
The definitions and properties of this section can be found in [TJ [22]. We begin with the 
definition of the anisotropic fractional Brownian sheet. 

Definition 1 (Fractional Brownian sheet) A fractional Brownian sheet (Wgl )< s ,t)E[0,l] 2 
with Hurst indices (a, (3) € (0, l) 2 is a centered two-parameter Gaussian process equal to 
zero on the set 

{(s,t) € [0,1] 2 , s = ort = 0}. 
For si,ti, S2,t2 £ [0, 1] the covariance function is given by 



R a ' f3 ((s 1 ,t 1 ),(s 2 ,t 2 )) := E 

= K a (s 1 ,s 2 )K"(t 1 ,t 2 ) 



-{s l +s 2 — |si — s 2 | J^(*i +t 2 -\ti-t 2 \ H ). 



We assume that (W s \ )( s ,t)e[o,i] 2 i s defined on a complete probability space (0, F, P) where 
F is generated by W a, P. Let us denote by T-L a 'P the canonical Hilbert space generated by the 
Gaussian process W a,/3 defined as the closure of the linear span generated by the indicator 
functions on [0, l] 2 with respect to the scalar product 

(![0, S i]x[0A]' 1 [0, S2 ]x[0,* 2 ])w«^ = ^^((Sl^l)) 0*2, *2))- 

The mapping l[o jS ]x[o,t] ^ W™t P rov ides an isometry between and the first chaos of 
W a, P denoted by . For an element cp of H a ^ we denote by W a '^((p) the image of ip in 
the space . 
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For any 7 in (0, 1) we denote by TV 1 the Hilbert space defined as the closure of the linear 
span generated by indicator functions on [0, 1] with respect to the scalar product 

(ilo.sib l[o, aa ])w = \ { s2 i + S V ~ l s i - s 2| 27 )- 

The space TV is in fact the canonical Hilbert space generated by the (one-dimensional) 
fractional Brownian motion with Hurst parameter 7 G (0, 1). With these notations we will 
often use the practical relation 

(l[0,si]x[0,ti], 1 [0,s 2 ]x[0,t 2 ])^.' 9 = (l-tQ.si], 1[0, S2 ])w q ( 1 [0,<i]' 1 [0,t 2 ]}#. V(si, S 2 , ti, t 2 ) G [0, l] 4 . 

More generally, for any two functions f,g G H 01 ^ such that Jj ^ 4 \f(u, v)g(a, b) \ \u — 
a |2a-2| t , _ b| 2 P~ 2 dod6dwdv < 00 we have 

(f,g) Ha .p =a(a)a{j3) [ f(u, v)g(a, b)\u - a\ 2a - 2 \v - b\ 2 ^- 2 dadbdudv (1) 

with a(a) = a(2a— 1). Note finally that we can also give a representation of (W^)(s,t)e[o,i] 2 
as a stochastic integral of kernels K a and if* with respect to a standard Brownian sheet 

(W( Sl t))(s,t)e[0,l} 2: 

W^ t) = J* ^K a (s,u)KP(t,v) dW M , (s,t) G [0, l] 2 , 

where K a is the usual kernel of the fractional Brownian motion B a which appears in its 
expression as a Wiener integral Bf = K a (t, s)dW s (see e.g. [18] for an explicit definition 
of this kernels; we will not use it in this paper). Using this representation, Tudor and Viens 
in [221 123j have developed a Malliavin calculus with respect to W a ^ . 

Let us recall the notion of self-similarity and stationary increments for a two-parameter 
process (see [2]). 

Definition 2 A two-parameter stochastic process {X 3) t)f 8 ,t)eT) T C R 2 has stationary in- 
crements if for every n G N and for every (si,t\), (s2j*2)j ■ • • ; (s n ,t n ) G T the law of the 
vector 

[X s+Sl>t+tl , X s+S2;t+t2 , . . . , X s+Sn;t+tn ) 
does not depends on (s,t) G T. 

Definition 3 A two-parameter stochastic process (X s j)( s ,t)eT> T C R 2 is self-similar with 
order (a, f$) if for any h, k > the process (X Syt ) s ,teT 

X st := h a k?X ± t 

' h'k 

has the same law as the process X. 

Note that the fractional Brownian sheet W a '@ is self-similar and has stationary increments 
in the sense of Definitions [3] and [2] (see [2]). 

Now we present some elements of Malliavin calculus with respect to fractional Brownian 
sheets and especially the Malliavin integration by parts formula ([4]). 
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2.2 Malliavin calculus for the fractional Brownian sheet 

We recall some definitions and properties of the Malliavin calculus for the fractional Brow- 
nian sheet. For general Gaussian processes these are contained in the framework described 
in [18]. 



By C£°(R ra ) we denote the space of infinitely differentiable functions from R n to R with 
bounded partial derivatives. For a cylindrical functional F of the form 

F = f(w°>( ) (<p 1 ),...,W a >f i (<p n j), n>l, t/?ij ... ,(p n E H a 'P, /eC 6 °°(R n ), (2) 

we define the Malliavin derivative DF of F as, 

n 

DF := £ dif (w^tpi), • • • , W a >P(vn)) Vi- 
i=i 

Furthermore D : L 2 (Q.,T , P) — > L 2 (Q,J 7 ,P;'H a 'P) is a closable operator and it can be 
extended to the closure of the Sobolev space D 1 ' 2 defined by the functionals F whose norm 
11-^111,2 is finite with 

\\F\\l 2 :=E[F 2 } +E[\\DF\\ 2 H ^] . 

The adjoint operator I\ of D is called the divergence operator and it is defined by the 
following duality relationship 

E[FI 1 (u)]=E[(DF,u) Ha ,,}, 

for F in B 1,2 and for u in T-L a, P such that there exists a constant c u > 0, satisfying 

\E[{DG ,u)y^a,fi]\ < c^HGH^a^ jrp), for every functional G of the form ([2]). 

Let n > 1. The nth Wiener chaos S) n of W a ^ is the closed linear subspace of L 2 (f2, T, P) 
generated by the random variables {H n {W a, P(tp)) , ip G T-L ' 13 , \\(p\\^a,p = 1} where H n de- 
notes the nth Hermite polynomial 

„ , , (-1)" si dl ( _*i 
HJx) = - — i—e 2 - — e 2 
v ; n! dxi V 

A linear isometry between the symmetric tensor product {H a, ^ Qn and S) n is defined as 

4 := n!if n (V a ^)) . (3) 

We conclude this section by the following integration by parts formula: 

E [FI n (h)\ = E 



(D n F,h) /n ,^^n 



where D"' 2 is the space of functionals F such that ||-F|| n ,2 is finite with 



he [H a A &n , F er 2 , (4) 



|F|| 2 j2 :=E[F 2 ]+^2E[\\D l F 



2 
i=l 
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2.3 Hermite variations of the fBs 

Let (W^l )s,t>o be a fractional Brownian sheet with Hurst parameter (a,/3) 6 (0, l) 2 . We 
will define the Hermite variations of order q > 1 of the fractional Brownian sheet by 

N-l AI-l 



iV — ± 1V1 — 1 / J 

Vn,m := E E ^ ( NaM ^ [W?£i±l- W fi±l - W< tl ± + W ± 

~Z ._ \ \ N < M N ' M N ' M N< 



i=0 j=0 

where H q is the Hermite polynomial of order q. Note that 



(5) 



W 



_ 2+1 
N ' M 



i+1 J- _L JL 



2a M -2/3 



N~ M 



which explains the appearance of the factor N a M 13 in ([5]): with this factor the random 
variable N a M? ( Wjf m - W^j+i ~ W%[ ± + W^j_ \ has L 2 -norm equal to 1. 

\ N > M N ' M N ' M N > M J 



We will use the notation 
i i + 1 



Ai 



N' N 



and Ai, j 



' i i + 1' 




" j 


i + i] 


_JV' N 


X 


_M' 


M 



Ai x Aj, 



for i G {0, . . . , N — 1}, j G {0, . . . , M - 1}. In principle Ai = depends on TV" but we 

will omit the superscript N to simplify the notation. With this notation we can write 



vv i+1 j + 1 



W a / J+1 - j + Wf P j = h ( lfi i+l lx rX i+i, j = /l(lAij) = /l(lAixAi) 

7T> M JV>~M~ TT'M iV'M V tN< N 1 L M ' M V 



r«,/3 



Here, and throughout the paper, indicates the multiple integral of order n > 1 with 
respect to the fractional Brownian sheet W a ^ . Since for any deterministic function h £ H a >P 
with norm one we have 

H q (h{h)) = -U^n 



we derive at 



N M , \ 

^ = ^EE™'% ^ixix ^ • 

5 1=1 i=l \ ljV> AT J x l M ' M V 

We want to study the limit of the (suitably normalized) sequence Vn,m as iV, M — > oo. 
Since this normalization is depending on the choice of a and /3, we will normalize it with a 
function (p(a, f3, N, M). 
Let us define 

N M , 

Vn,m := -<f(a, (3, iV,M) EE M X l 
q - i=i j=i ^ 

By renormalization of the sequence V/v,m we understand a function if (a, j3, N, M) to fulfill 



_L i+Ii v r J- i+1' 

L AT ' AT J L M ' M 



(6) 



the property KVj^ M N,A £^ 00 \ 
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It turns out that the limit of the sequence Vn m is either Gaussian, or a Hermite random 
variable, which is the value at time (1, 1) of a two-parameter Hermite process. 
In the case when V/v,Af converges to a Gaussian random variable, our proof will be based 
on the following result (see [15], see also [IS]). 

Theorem 1 Let F be a random variable in the qth Wiener chaos. Then 



d(F,N) < cJE((l-q- 1 \\DF\\l ChP ) (7) 



The above inequality still holds true for several distances (Kolmogorov, Wasserstein, total 
variation or Fortet-Mourier) . The constant c is equal to 1 in the case of the Kolmogorov 
and of the Wasserstein distance, c=2 for the total variation distance and c = 4 in the case 
of the Fortet-Mourier distance. 

Until the end of this paper d will denote one of the distances mentioned in the previous 
theorem. We will also assume that q > 2 because for q = 1 we have Hi = x and then 
Vjv,m is Gaussian; this case is trivial. Our argumentation has the following structure. We 
first compute the Malliavin derivative (with respect to the fractional Brownian sheet W a, P) 
DVn,m and we compute its norm in the space Ti a ^ . We will get 

1 N M , v 

DV N>M = 7TT¥>(a!, P, N,M)^2^2 Iq-1 iff il lx[ i i+J]x[i i±il, 

{q-l)l , , V Wm'TvtJ/ [ n> n I*Im' m 1 

i=l j=l x ' 

and 

. AT-1 M-l 

\\DV NjM \\l_ a ,, = _ (<p(a,p,N,M)f x Yl ( 1 Aw(-).Wj'('))«^ 

i,i'=0ij'=0 

The product formula for multiple integrals (see [18], Chapter 1) reads 

Vi (l K) = EpI(C^ 1 ) a <lAy(-),WjK-%^ 

p=0 



Xlo 



Sg-1-pA.-l ®g— 1-p 
2g-2-2p I J-Aij ® A Ai',.j' 



where Cg_i : = ^ ^ for q > 2,p < q — 1 and /<§>5 denotes the symmetrization of the 
function f ® g. Hence, we have 

j iV-l M-l 9-1 

IIW^a/II^ = _ (y(a,/3,iV,M)) 2 x ^ £ X>AijO,WjK0>K>pl 

^ i,i'=0 jj'=0p=0 

X (Cq_i) 2 hq-2-2p (l/ 



7— 1— 3»,Svi ®?-l-p 
L Aij Q9 - L Ai',j' 
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Let us isolate the term p = q — 1 in the above expression. In this case 2q — 2 — 2p = and 
this term gives the expectation of H-DV^m!!^*^- 

||-DVjV,Af||^a,/3 (8) 

' : ( V (a,f3,N,M)) 2 



(q-l)l 2 

N-l AI-l q-2 

x E E E( 1 A^(0,iA^iK0)SV!(^_j 2 /2 ? -2- 2 p(iif- 1 -^ii^- 

i,i'=0 j,j'=0 p=0 

N-l M-l 

+ T—jy^,P,N,M)) 2 ]T £ (l Ai)i (.) ; W,i'(-)>^,. =:T 1+ r 2 . (9) 

i,i'=0j'j'=0 

The term T 2 is a deterministic term which is equal to E||DVa/" ) m||Lc«,/3- 
With the correct choice of the normalization we will show that T 2 is converging to q as 
N, M goes to infinity and T\ converges to zero in 1? sense. Using Theorem [1] we will prove 
the convergence to a standard normal random variable of Vn,m an d we give bounds for 
the speed of convergence. The distinction between the two cases (when the limit is normal 
and when the limit is non-Gaussian) will be made by the term T\: it converges to zero if 
a < 1 — 7^or/3<l — 7^, while for a, (3 > 1 — ^ this term converges to a constant. 
Let us first discuss the normalization <p(a, (3, N, M) and the convergence of T 2 in the fol- 
lowing lemma. Given two sequences of real numbers (a n ) n >i and (b n ) n >i, we write a n < b n 
for sup^ j^f < oo. 

Lemma 1 Let T 2 be as in ([9]). Then q~ 1 T 2 N,R Ij^ 00 \ f or ^ e following choices of ip: 

1) (p(a t p,N,M) = J^N aq ~ l l 2 M aq - l l 2 , ifO< a,/3 < 1-i and q" 1 ^ - 1 < N' 1 + 

jy2qa-2q+l _|_ j^-1 _|_ ^2g/3-2g+l ^ 



2j <p(a,/3,N,M) = ^^A^-i M9 -i( logM) -i/2 ; z/0 < Q < j _ ^ ^ = x _ i_ and 
qT X T 2 - 1 < N- 1 + N 2qa ~ 2q+1 + (logM)- 1 , 

5; tp(a, (3, N, M) = N q ~ x (log Ny^M*- 1 (log M)- 1 / 2 , if a = (3 = I - ± and 

q-±T 2 - 1 < (logA^)- 1 + (logM)- 1 , 



4) <p(a,(3,N,M) = J^L-N^-WM*- 1 , if < a < I - ±, (3 > I - ± andq- l T 2 -l< 

N -l + N -2qa+2q-l + M ~2q(5+2q-l 
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5) <p(a,P,N,M) = N*- 1 (log Ny^M*- 1 , if < a = 1 - ^, > 1 - £ and 
qr x T 2 - 1 < (log(A^)- 1 + M~ 2qfS+2q ^ 1 . 

6) ip(a,P,N,M) = ^Zm- l Mi- 1 , if a > 1 - i /3 > 1 - i and q~ l T 2 - 1 < 

N -2qa+2q-l + M -2q/3+2q-l _ 

where s., i. and k. are defined in Lemma\^ in the Appendix. 

Proof: Using the properties of the scalar product in Hilbert spaces we have 
. iV-l Af-l 

T2 = 7— w^'^^^) 2 x E E^wo^wO)^ 

i,i'=0j,j'=0 



/JV-1 \ /M-1 

-— M«,/3,iV,M)) 2 x ^(i A i(0,W(-)>V x E (WO.WO) 
^' \i,i'=o / \iJ'=o 
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The result follows then from Lemma [2] in the Appendix. I 

Remark 1 As mentioned above, T 2 = K\\DVnm \\h<*,p- On the other hand, we also have 

qT 2 = EV 2 >M . 

Indeed, this is true because for every multiple integral F = I q (f), it holds that EF 2 = 
qE\\DF\\ 2 Ha ^. 

3 The Central Limit Case 

We will prove that for every a, 13 G (0, l) 2 \ ( 1 - 2 1 , l) a Central Limit Theor em holds 



2g> 

where Vn,m was defined in ([6]). Using the Stein's method we also give the Berry-Esseen 
bounds for this convergence. 

Theorem 2 (Central Limits) 

Let Vn,m be defined by |]5p. For every (a, ft) G (0, l) 2 , we denote by c a ^ a generic positive 
constant which depends on ct,/3,q and on the distance d and which is independent of N and 
M. We have 

1) If < a, j3 < 1 — then Vn,m converges in law to a standard normal r.v. M with 
normalization <p(a,f3,N,M) = ^ j^N aq ~ 1 / 2 M aq ~ l l 2 . In addition 

d{V NM M) < Ca^VN- 1 + N 2a - 2 + N 2aq - 2q + l + M~ l + M 2 ?- 2 + M 2 P q ~ 2q+l . 
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2) 7/0 < a < 1-^ and (3 = 1-^, thenV N . M converges in law to a standard normal r.t 



N with normalization <p(a,(3,N,M) = J ' -^N aq ~ l W 1 * 1 (log M) -1 / 2 . j n addition 
d{V NM M) < c a ,^N~ l + N 2a ~ 2 + N 2 "*- 2 ^ 1 + (log M) _1 . 

Af converges in law to a standard normal r.v. J\f with 



normalization ip(a,(3,N,M) = J^-iV 9-1 (log N)~ 1 / 2 M q - 1 (log M)" 1 / 2 . 7 n addition 
d(V NjM ,Af) < c^VlogA^" 1 +logM" 1 . 

If a < 1 — ^- and j3 > 1 — i/ien V/v^f converges in law to a standard normal r.v. 
N with normalization (p(a,/3,N,M) = J ' j^N aq - 1 / 2 M q ~ x . In addition 



d{V NM M) < C^VN- 1 + iV 2 "" 2 + jV2/3<7-2,+ l + M 2/3«,-2,+ l. 

5j 7/0 = 1 — 5- and /3 > 1 — 5-, i/ien Vjv^ converges in law to a standard normal r.v. 
N with normalization (p(a,/3,N,M) = ^J T ^N q ~ 1 (log N)~ 1 / 2 M q ~ 1 . In addition 

d(V NM ,M) < c a ^log(N)-^ + M 2 ? q - 2q +\ 

Proof: Recall that 

||7JVfyAf = -Ti + T 2 , 

where the summands T\ and T 2 are given as in ([9]). We apply Lemma Q] to see that 1 — q~ 1 T 2 
converges to zero as N, M goes to infinity. 

Let us show that 7\ is converging to zero in L 2 (S7). We use the orthogonality of the iterated 
integrals to compute 



(9 - 

yH ' i,i',k,k'=0j,j',t,e'=0p=0 
( C 'g-l) 4 ( 1 Aij(-)) 1 Ai',i'(-))^^( 1 AM(')) iAfe'^KO)^,^ 



XE 



hq-2-2p (l Ai j ^Ai'J' ) 7 2g-2-2p ( 1 Afc,^ (g)1 Afc'/' 



JV-1 Af-1 g-2 
? — -^(a,/3,iV,M)) 4 £ £ £(p!) 2 
w i,i',k,k'=oj,j',e,e'=op=o 

(Cg_ 1 ) 4 (lAij(0> lAi' ) j'(-)}^^( 1 AM(") ) lAJfeV'(O)^^ 
vf9« 9 9r.M/l lS>9 ~ 1 ~ P (9i1 <g ' 9 ~ 1 ~ p 1 ^S-l-P^-i 8>?-l-P\ 
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Now, let us discuss the tensorized terms. We use the fact that 
i ®g-i-p<>-i ®5-i-p i ®?-i-p,Svi ®?-i-p 



a+6=g— 1— p;c+d=q— 1— p 

X (lAi'J', lAf;i}^,(3 (lAi'j', lAfc'^'}^a,^ 

(lAt, lAfc)H a ( :1 -Aj, 1a^)^/3 

a+6=g— 1— p;c+d=q— 1— p 

x(1Aij lAfc'^^Ajf, lA^')^(lAi'j 1a*)k«(1Aj'j 1a^)^/9 



x (lAi' 5 lAfe')'H a ( : '-Ai', lA£'}-H/3 



(we recall that 1a« : = ljj. i±ii)- Therefore we finally have 

1 

(^1) 



p=0 a+6=g— 1— p;c+d=g— 1— p 

a N (p, a, a, b, c, d)b M (p, P, a, 6, c, d), 



with 



iV-1 

a,N(p,a,a,b,c,d) = (lAi, lAfc)ft«(lAi, lAfc')tt<* 

i,i',k,k'=0 

(lAi'j lAfc)^a(lAt') lAfe')w a 
(lAi, lAi')?£« (lAfcj lAfc')^t 

and 6m(P) /3,a,b,c,d) similarly defined. We apply Lemma [3] stated in the Appendix to the 
terms a^v and 6m to conclude the convergence of Ti to zero. Hence, E[(l— g^ll-DV/v.Afll^^) 2 ] 
<7 _2 E[|T 1 | 2 ] + (1— g~ 1 T 2 ) 2 which converges to zero for a < 1 — ^or/3 < 1— The bounds on 
the rate of convergence are given by the Lemmas [3] and [TJ Using Theoremfll the conclusion 
of the theorem follows. ■ 

The fact that the term T\ converges to zero makes the difference between the sit- 
uations treated in the above theorem and the non-central limit case proved in the next 
section. 

4 The Non Central Limit Theorem 

We will assume throughout this section that the Hurst parameters a, P satisfy 

1 > a,P > 1- — . 

2q 
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We will study the limit of the sequence Vjv m given by the formula ([6]) with the renormal- 
ization factor ip from Lemma [TJ point 6. Let us denote by /tjv,M the kernel of the random 
variable Vn m which is an element of the qth. Wiener chaos, i.e. 



JV-1M-1 

h N , M = ^(a,/3,iV,M)^ 

i=0 ]=0 

We will prove that {hpfj{)N,M>l is a Cauchy sequence in the Hilbert space (n^f q . Using 
relation (JTJ), we obtain 

N-lN'-i / ,i±i ,£±i \ q 

u-v 2a - 2 dudv 



(a(2a - !))« E E /, " /„" 

i=0 i'=0 \ 7f W 
M-lM'-l ( r i±l 

W-i))'E E // // k-h 2/3 - 2 ^ 
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j=0 j' = l \ M M 

and this converges to (see also [5] or 



C2(a,/3)^(a(2a-l)) 9 (/3(2/3-l))« / f \u - v^^Hudv f f \u - v^'^dudv, 
1- Jo Jo Jo Jo 

where C2(a, 8) = The above constant is equal to 

c 2 (a,f3)^-(a(2a - l)) 9 (/3(2/3 - l)) q - -— ^— 

1 ,HJ ql 2K V W l ^ ^ " (aq-q + l)(2aq-2q+l) (Pq-q + l)(2Pq-2q + l) 



Remark 2 Note that the above constant is actually . 

It follows that the sequence Tijv.Af is Cauchy in the Hilbert space (?i a 'P)® q and as N, M — > oo 
it has a limit in (l~L a, P^ q denoted by pS- q \ In the same way, the sequence 

[(iV-l)t] [(M-l)s] 

h N , M (t,s) = - V (a,P,N,M) E E ^iiiiixM-itil 

y! — „ — " ljv> jv I a Im' m ) 
i=0 j=0 

is Cauchy in {ji a >P^® q for every fixed s,t and it has a limit in this Hilbert space which 

will be denoted by fj, q j. Notice that (1,1) and that is a normalized uniform 

measure on the set ([0, t] x [0, s]) 9 . 
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Definition 4 We define the Hermite sheet process of order q and with Hurst parameters 
a,(3 G (0,1), denoted by {Zffl^o^, by 

zlf^:=Z$=I q ($), VM€[0,1]. 
The previous computations lead to the following theorem. 

Theorem 3 Let Vn,m be given by ^ with the function <p defined in Lemma [JJ, point 6. 
Consider the Hermite sheet introduced in Definition^ Then for q > 2 it holds 



lim E[\V NM -Z\ 2 } = 0, 

N,M->°o 



where Z := Z^\. 



Proof: Note that ^E[\ V NiM -Z\ 2 ] = ||^Ar,M||( Ha ,/3 ) ® ? + ||^ (9) ll( W c l ^)^-2(/iAr>/, / u (9) )( WQ , / 3)(g, g . 
The computations of the beginning of this section complete the proof. ■ 

Let us prove below some basic properties of the Hermite sheet. 

Proposition 1 Let us consider the Hermite sheet (-^i t ) S| tg[Q,i] from Definition^ We have 
the following: 

a) The covariance of the Hermite sheet is given by 

— Rq{a-l)+l( s , u)R q (/3-l)+l{t,v). 

Consequently, it has the same covariance as the fractional Brownian sheet with Hurst 
parameters q(a — 1) + 1 and q((3 — 1) + 1. 

b) The Hermite process is self-similar in the following sense: for every c,d>0, the process 

z { $ •= (z^U dt 

has the same law as c q ( a -^ +l d q ^~^ +1 Z*fJ . 

c ) The Hermite process has stationary increments in the sense of Definition [H 

d) The paths are Holder continuous of order (a',/3') with < a' < a and < /3' < j3. 
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Proof: Let / be an arbitrary function in {H a ^^ q . It holds that 



" I- AT' JV l X Ul' M 1 \ / 



W-^Mi- 1 [(N ~ 1)t] [{M ~ 1)S] 



a(a)ia((jyc 2 (a,f3y 1 / 2 { £ £ 



/ dxi..dx q dyi..dy q f ((xi,yi), ..,(x q ,y q )) x / dai..da g dbi..db g 



'[0,1] 29 

xnia,-x,r- 2 ni6 fc - yfc r- 2 



— > a 



fc=l fc=l 

■t 



(a) q a(j3) q C2(a, ft) 1//2 — - / cia / rft> / dx\..dx q dy\..dy, 



Q ] - Jo Jo J[o,i] 2 i 



c/((xi,yi),..,(x g ,y 9 )) Yl \a-x k \ 2a 2 \b-y k f 



1 2/3-2 

k=l k=l 



By applying the above formula for / = ^1 and using the fact that 



we obtain the point a). 
Concerning b), let us denote by 



[(JV-l)t] [(M-l)]s 

H N>M (t,s) = - ] <p(a,f3,N,M) £ £ I g (iff ^ m ) 

i=0 j=0 



We know that 

H cN4M (t,s) '-► Zfl (10) 



in L 2 (0) for every s,t € [0, 1]. But 

[I.V-J)/. (.U-Ji 

(ciV) 1 -( 1 - Q )9(dM) 1 -( 1 -/ 3 )'? 



i=Q j=0 



H NM (t,s) -> Zf> (11) 



( c )l-(l-a) 9( - d )l-(l-^g 



The point b) follows easily from (|10|) and (|lip . 

Point c) is a consequence of the fact that the fractional Brownian sheet has stationary incre- 
ments in the sense of Definition [2] while point d) can be easily proved by using Kolmogorov 
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continuity criterion together with points b) and c) above (see also Section 4, page 35-36 in 
[2j). ■ 

5 Appendix 

We recall the following two technical lemmas which have been proved in |15| and [5]. 
Lemma 2 Let 7 in (0, 1) and q be an integer with q > 2. We set 

r-y(z) := X - (\z + 1| 2 ^ + \z- 1| 2 ^ - 2|z| 2 ^) , z € Z. 

We have: 
(i) IfO < 7 < 1 - then 



2q 

V- I 

JV->oo 



lim N 2 ^- 1 (lAi, =E r 7^) 9 = : s 7> 



i,i'=0 reZ 
and liV 2 ^" 1 ESiaAi, lAi')|f7 " *yl < N' 1 + iV 2 ^" 2 ^ 1 . 
^ 7/7 = 1--^, t/ien 

jim log^)- 1 ^" 2 £ (l Ai , W)^ = 2 A2g- l)(g-l) V =: 
i,i'=0 \ 9 / 

and I log^)" 1 ^" 2 EffioilAi, W>^ " ^| < log(iV)- 1 . 
(m,) 1/7 > 1 — j^, t/ien 

iV-1 



^ U { lq -q+l)(2 lq -2q+l) 



and -a^I^JV^- 1 - 2 ^. 

Proof: The first two claims can be found respectively in [T5J p. 102], O p. 491-492]. For 
the third part we define f N := N^ 1 lf[ m Then is a Cauchy sequence in 

(W)® q with limit / and ||/|| 2 # 7 )® 9 = K-y. For the rate of convergence we have 

II/^II^t)® 1 ? ~~ II/II(H7)®9 = II/JV — /l| 2 fl-7)®g + 2(/jV — /, f)(Hf)®9 

< ||/jV ~~ /Il(jf7)®9 + 2||/jV — /|| 2 JJT)®9 H/ll^iTr)®?" 
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Refer to [5j Proposition 3.1] to see the details and to get that the order is 0(N 2q 1 2jq ) (a 
direct argument as in the proof of the next lemma can be also employed). 

■ We also state the following estimates which have been obtained respectively in [151 P- 102, 
p. 104] and in [5, p. 491-492]. 

Lemma 3 Let 7 in (0, 1). We set q,p, a, b, c, d integers such that: q > 2, p € {0, . . . , q — 2} 
and a + b = c + d = q — 1— p. We have: 

(i) If < 7 < 1 - ^, then 

N-l 

_/y4 9 7-2 ^2 (Iai, lAi')^ (lAfc, lAA')?£7 (IAi, 1 Ak)u~t {1 Ai' , ^Ak')nt 

i,i',k,k'=0 

x(lAi') lAfc)^7(lAi' 5 1~Ak')u~t 
< AT" 1 + iV 2 7-2 + ^y2 7 g-2g+l^ 

fWy) 7/7 = 1--^, i/ien 



1 f_/\T) 2 ^ (^Ai, lAi0^ 1 (lAfc, lAfc')?^ 1 ( :1 -Aj, lAfc)^( 1 Ai' ; lAfc')«7 
°^ ' i,i',k,k'=0 

X (lAi'j lAfc)«T (lAi', lAfc')w 
< log(Af)^ 1 . 

(raj If j > 1 — then 

N-l 

jyiq-A (lAi, lAiO?^ 1 (lAA:, lAfc')^ 1 (lAi, lAfc)?{7 (lAi', lAfc'}«7 

i,i',jfe,Jfe'=0 

x (lAi', lAfc)^7(lAi', lAfc')^7 

< 1. 



Proof: The first point is proved in P2J p. 102, 105]. The point (ii) is done in p. 491-492]. 
The last case can be treated in the following way. The quantity (Iaij ^-Ai')w is equivalent 
with a constant times N~ 2ry \i — i'\ 2 ^~ 2 and the sum appearing in (iii) is then equivalent to 

N-l 

N 4q-4 N -4 lq \{ _i'|( 2 7-2)(p+l)| fc _ fe / 1 (2 7 -2) (p+1) 

i,i',k,k'=0 

x \i - k\^~ 2 ^\i' - k'\^- 2 >\i' - k\^- 2 ^\i - k'\^- 2 ^ d 
N-l 

= N ~* N- 2q ^- 2 ^\i-i'\ ( - 2 ^- 2 ^ p+ ^\k-k'\ ( - 2 ^- 2 ^ p+ ^ 
i,i',k,k'=0 

X |- _ A .|(2 7 -2)a|./ _ ^[(27-2)6^/ _ fc|(2 7 -2)c|- _ k '\(2y-2)d^ 
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for N large enough and this is a Riemann sum which converges to a constant. 
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